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Mineralized biological materials such as bone, sea sponges or diatoms provide load-bearing and armor 
functions and universally feature structural hierarchies from nano to macro. Here we report a systematic 
investigation of the effect of hierarchical structures on toughness and defect- tolerance based on a single and 
mechanically inferior brittle base material, silica, using a bottom-up approach rooted in atomistic modeling. 
Our analysis reveals drastic changes in the material crack-propagation resistance (R-curve) solely due to the 
introduction of hierarchical structures that also result in a vastly increased toughness and defect-tolerance, 
enabling stable crack propagation over an extensive range of crack sizes. Over a range of up to four hierarchy 
levels, we find an exponential increase in the defect-tolerance approaching hundred micrometers without 
introducing additional mechanisms or materials. This presents a significant departure from the 
defect-tolerance of the base material, silica, which is brittle and highly sensitive even to extremely small 
nanometer-scale defects. 



This work addresses a fundamental question — is it possible, and if yes, how, to create materials with enhanced 
mechanical properties based on a single and mechanically inferior material, solely by using geometry as a 
design parameter? Many biological materials with a structural purpose, such as bone, nacre, sea sponge 
exoskeletons, diatoms and spider silk' showcase the use of structurally inferior building blocks, arranged in 
multiple levels of hierarchical structure. These materials are used for load-bearing, skeletal support or armor 
protection applications and require excellent mechanical properties, including great robustness against cata- 
strophic material failure. In particular, high toughness-quantifying the capacity of a material to dissipate energy 
under large loading-is a vital material property for these biological applications. Experimentally, the fracture 
behavior of materials such as bone and nacre have been shown to be intimately linked across multiple structural 
length scales'"'", with mechanisms at several different hierarchy levels participating in the overall behavior'^ A 
fundamental question, however, remains, whether or not specific base materials or mechanisms (e.g.: reinforce- 
ment materials, an organic phase, specific interfaces or interfacial effects, grain boundaries, etc.) are needed in 
order to achieve enhanced mechanical properties at larger scales. Since the impact of hierarchical structures alone 
on toughness cannot easily be quantitatively compared to other toughness improvement techniques, the under- 
standing and use of hierarchies as a sole design variable in the bottom-up design of materials remains limited. 

A thorough understanding of the mechanical properties of hierarchical materials can be achieved through the 
development of models that appropriately reflect their structural and mechanistic nature from the bottom up. A 
pioneering theoretical approach in addressing this issue has been the study of self-similar hierarchical assemblies 
with some assumptions about similarity in failure mechanisms at each hierarchy level, using continuum 
mechanics'^"*. These continuum-level models can predict the strength, stiffness and toughness scaling with 
the number of hierarchy levels and have been applied to several cases. Applications included a self-similar 
hierarchical assembly of the bone nanostructure"" and gecko adhesion'"'™. Self- similarity in geometric design, 
however, is not always observed in biological materials, and in particular in bone, nacre, sea sponge exoskeletons, 
diatom algae and other case^'*'^''^^, the structural arrangement at each hierarchy level is quite different. Another 
approach has been to take into account explicit probabilistic modes of failure pathways in several levels of 
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hierarchy, under the assumption of certain unit failure events (e.g., 
failure of unit hydrogen bond clusters in hierarchical alpha-helfac 
structures)^'''^'*. Some qualitative models have been derived to identify 
design similarities in the hierarchical structures of biological 
materials that may offer the key to improvements of properties^^"^^. 
Limitations to existing models remain with respect to their ability to 
provide a complete description of the mechanism of failure for 
example via crack propagation. 

Here we use a bottom-up approach to systematically examine the 
effects of hierarchical structures on fracture toughness by investi- 
gating how crack propagation is modified solely through the pres- 
ence of hierarchical structures and without adding additional 
materials or mechanisms. We specifically consider that, from a frac- 
ture mechanics^" point of view, a single major crack propagates 
through the material and leads to failure and that therefore this one 
crack interacts with all levels of hierarchy. We also examine the 
defect-tolerance length scale, which measures the sensitivity of 
fracture strength to the presence and size of a crack; where a higher 
defect-tolerance implies a lower sensitivity to crack size, i.e., a large 
change in the size of a crack in the material leads only to a relatively 
small change in the fracture strength. 



We use silica as a single base material for constructing various 
hierarchical structures to represent a set of model materials that 
are used to carry out a series of computational experiments. The 
choice of silica as a base material is inspired by the dominant material 
found in many mineralized organisms such as diatoms and deep sea 
sponges'^'^', which showcase the use of a predominantly sUica-based 
structure (with up to 97-99% silica content) that result in protective 
exoskeletons that perform much better than the base material alone; 
which is stiff but very fragile and thus mechanically inferior. 

Previous work based on first principles based reactive modeling 
has demonstrated that nanoporous silica as found in diatom algae, 
for example, is mechanically deformable and "ductile" (tough), but 
very soft'^''''' (Fig. la-b). It was demonstrated^"* that the geometry of 
this structure, with nanoscale silica struts confined to a few nan- 
ometer cross-section, is critical in providing the highly deformable 
material behavior; which also agrees qualitatively with size-depend- 
ent experimental results in silicon'^. The earlier atomistic modeling 
showed that each strut is in a flaw-tolerant state that is characterized 
by a largely homogeneous stress state", facilitating a ductile-like 
behavior of silica. Even though the individual struts reach very high 
strengths under this geometric confinement, the overall strength of 
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Figure 1 Structure and properties of nanoporous silica inspired from the nanostructure of diatom algae, and two commonly occurring periodic 
motifs chosen for the nanoporous silica/bulk silica composites, a. Scanning Electron Microscope (SEM) image of a diatom nanostructure (image 
reprinted with permission from Ref. 49) and fuU-atomistic model of a-quartz nanoporous silica, stress-strain response obtained from atomistic 
simulations of tensile loading of the nanoporous silica and bulk silica, b, Mesoscale model stress-strain response. The base material for both structures is 
silica and identical for both structures but the different organization leads to a drastically changed mechanical response (the stress-strain data shown is 
averaged over a representative cell in each case) . c. Representation of two biological structures, hone-like and biocalcite-like. The upper structure (hone- 
like) represents one in which the soft/tough nanoporous silica material is the matrix phase in which hard/brittle bulk silica platelets are dispersed. In the 
lower structure (biocalcite-like) the hard/brittle bulk silica material serves as the matrix material in which soft/tough nanoporous silica platelets are 
dispersed. These structures mimic those seen in bone and nacre biological calcite single crystals, respectively'^""*. 
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the nanoporous material (averaged over a representative element of 
the porous material) is smaller due to the relatively low volume 
fraction of silica (Fig. lb). 

The natural nanoporous geometry found at the lowest level of the 
structural hierarchy, and the corresponding nanoporous a-quartz 
silica structure inspired from this geometry, are shown in Fig. la. 
The soft, deformable and "ductile" (tough) nanoporous silica is com- 
bined with stiff and brittle bulk silica, forming two material phases as 
fundamental building blocks at the next level of hierarchy that allows 
us to use a single material as the fundamental building block, but 
organized hierarchically in different geometries. Many biological 
materials contain not only mineral phases but a combination of 
mineral and organic (e.g.: protein) phase. In the spirit of using a 
single material in constructing hierarchical materials, the protein 
phase found in biological materials such as bone, which is soft, 
extensible and "ductile" (tough), is replaced by nanoporous silica 
in our structures. The mineral constituent, which is hard and 
brittle, is replaced by bulk silica. No interfacial effects between any 
two phases are considered, reflecting that all structures are 
composed of a perfectly bonded single material as the fundamental 
building block. 

To enable us to access the scales associated with such hierarchical 
materials (approaching hundred micrometers) we use an in silico 
multiscale approach based on a mesoscopic model that describes 
the material as a collection of particles connected by nonlinear 
springs (details of the atomisticaUy-informed mesoscale method 
are presented in the Methods section, the Supplementary 
Information material, as well as in'"'). We emphasize that the model 



used here is designed by our desire to derive generic insight into the 
mechanical properties and deformation mechanisms valid for a 
broader class of biological materials, and not for one specific 
material. While such a simple model formulation does not allow us 
to derive quantitative properties of specific types of biological mate- 
rials, it enables us to understand generic mechanisms that control the 
failure behavior of hierarchical materials and to carry out a system- 
atic comparison of different structural designs. 

Results 

Two-hierarchy level structures. We first consider regularly 
distributed composite structures with two levels of hierarchy. For a 
particle-reinforced composite, two distinct design schemes are 
possible. One in which the soft material is the matrix, and another 
one where the soft material makes up the reinforcement particles. 
Both cases are observed in biological materials, and in order to 
investigate both we choose two representative systems. In the first 
one {"bone-like") the soft, deformable and "ductile" (tough) material 
is the continuous phase, with the hard and brittle material dispersed 
in platelet form within the matrix, mimicking the periodic 
arrangement observed in several hierarchy levels in bone and 
nacre'. The second case {"biocalcite-like") is where the hard and 
brittle constituent is the continuous phase, with regions of soft, 
deformable and "ductile" (tough) material embedded within the 
hard matrix. These structures are found in biological crystals 
where protein material is encapsulated in a hard crystal, such as 
biological calcite single crystals"'""*. The use of a periodic 
distribution of one material in another not only allows us to 




0.005 0.01 
strain 



0.005 0.01 
Strain 



0.015 



1 


1 






1 












1 


1 






Figure 2 | Characteristic stress-strain curves with and without a crack for two geometries and associated crack propagation paths (initial crack: white 
line, crack path: red line), a. Bone-like composite structures, with and without presence of a pre-crack (pre-crack length 5.4 ^m). The fracture strength 
changes drastically upon the introduction of a crack, b. Stress-strain curves for a biocalcite-like composite structure, with and without presence of a 
pre-crack (pre-crack length 6.96 [im). The sensitivity to fracture strength is much smaller than for the bone-like composite, although the magnitude of 
the fracture strength is lower, c and d, Crack pathways (marked in red) for hone-like and biocalcite-like hierarchical structures in the presence of a pre- 
crack. c. The image shows that for bone-like structures, the pre-crack propagates through the sample, but the structure is toughened by the platelets 
bridging the wake of the crack as it propagates, d. For small crack sizes failure in the biocalcite-like structure propagates through the nucleation of several 
micro-cracks at nanoporous silica/bulk-silica interfaces located far from the original crack tip. The fracture strength is reached when several of these 
micro-cracks link up along with the pre-crack to create a complete failure path through the sample. These results show that different crack propagation 
paths in the two structures lead to different defect-tolerance responses. 
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capture actual design morphologies observed in biological structures, 
but also removes stochastic elements of the dependence of 
mechanical properties on random morphologies. The volume 
fraction of the stiff silica phase is kept at a high value and constant 
at 80% for all structures under consideration; where this high value is 
important to maintain a high modulus. 

Both composite structures considered in the analysis of two- 
hierarchy structures are shown in Fig. Ic. For the bone-like arrange- 
ment, the platelets of bulk silica have a rectangular shaped 
cross-section and 8.4 |im by 2.4 |im in size in the X-Y plane, and 
100 nm in the out-of-plane direction. They are arranged in a stag- 
gered fashion, with an overlap equal to half their length across sub- 
sequent layers. The overall size of the structure considered here is 
27 [im by 70 |im and 100 nm out-of-plane. For the biocalcite-like 
arrangement, the soft nanoporous silica is embedded as rectangular 
inclusions with a cross-section size of 8.7 \xm by 0.7 |im in size in the 
X-Y plane, and 100 nm in the out-of-plane Z-direction. Next we 
create sharp edge cracks in all structures and load them under 
quasi-static tensile mode I loading. Initial crack sizes considered 
range from =5 to 20 |j,m, and crack initiation is identified by the 
advance of the crack front at a particular applied strain. Figure 2a 
shows representative stress-strain responses for structures with bone- 
like and biocalcite-like morphologies, with and without the presence 
of the edge crack. The data shows that both structures feature a drop 
in fracture stress under the presence of a crack. The bone-like struc- 
ture has a higher fracture stress but shows a larger drop in strength 
when a crack is introduced, whereas the biocalcite-like structure has 
smaller fracture strength, but also a lower sensitivity to the presence 
of a crack. Thus the biocalcite-like structure shows a higher defect- 
tolerance at the cost of a lower fracture strength. 
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Taking a closer look at the fracture mechanics and crack propaga- 
tion behavior in both structures. Fig. 2b provides an indication as to 
why the biocalcite-like structure feature a higher defect- tolerance. For 
small crack sizes, failure in these materials propagates through the 
nucleation of several micro-cracks at nano-porous sUica/bulk-silica 
interfaces located far from the original crack tip. This makes the 
failure morphology appear quite similar in the presence of smaller 
cracks and in the absence of any, and provides comparable fracture 
strength values. On the other hand, for the bone-like structures, the 
pre-crack always propagates as an individual crack, but toughening 
occurs by platelets bridging the wake of the crack as it propagates. 
This analysis shows that a change of the structural arrangements of 
the material has the potential to direct a significant alteration of the 
fracture properties, despite the same single material constituent and 
identical volume fraction in both cases. 

Extension to three- and four-hierarchy level structures. To assess 
the effect of additional hierarchy levels on the mechanical properties, 
the two-hierarchy structures discussed in the previous section are 
now extended to feature additional levels of hierarchy. Both self- 
similar geometries, in the vein of earlier studies"", and dissimilar 
geometries are considered (see insets in Fig. 3). Even though self- 
similar (fractal) geometries are not found widely in biological 
materials, they have been used in the literature to build theoretical 
and computational models of hierarchical structures and are 
therefore also included here. Dissimilar geometries are considered 
as representative cases of how the geometry at different hierarchy 
levels in bone, diatoms and other biological materials can be quite 
different from one another. The self-similar geometry uses a replica 
of the bone-like arrangement at different scales (Fig. 3b inset). The 
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Figure 3 | Comparison of the stress-strain response of two-hierarchy (a, b) and three-hierarchy (c, d) materials, considering both self-similar (a, c) and 
dissimilar (c, d) designs, a. Stress-strain curves for the two-hierarchy bone-like composite structures, with and without presence of a pre-crack (pre-crack 
length 5.8 nm). b. Stress-strain curves for two-hierarchy biocalcite-like composite structures, with and without presence of a pre-crack (pre-crack length 
6.96 |J.m). c. Three-hierarchy self-similar structure made of hone-like composite structure at both the second and third levels, with and without presence of 
a pre-crack (pre-crack length 4.8 pm). The data shows that the sensitivity of fracture strength vs. crack size is much smaller for the three-hierarchy 
material, d. Three-hierarchy dissimilar structure made of biocalcite-like composite structure at the second level and bone-like at the third level, with and 
without presence of a pre-crack (pre-crack length 5.8 ^m). Similar as for the case presented in panel c, the sensitivity of fracture strength vs. crack size is 
smaller for the three-hierarchy material than for the two-hierarchy case. 
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platelets size at the second level of hierarchy is 5.95 |im by 1.16 |im in 
cross-section, while at the third level it is 12.7 |im by 5.4 |im. The 
overall sample size is 54 |xm by 70 |j,m in the X-Y plane, and 100 nm 
in the out-of-plane Z-direction. The dissimilar geometry is based 
on the bio-calcite template for the second level, and bone-like 
arrangement at the third level of hierarchy (Fig. 3d inset). Here, 
the platelets size at the second level of hierarchy is 2.02 |.im by 
0.47 |j,m in cross-section, while at the third level is 12.6 |.im by 
5.5 ^m. Initial crack sizes in these models range from =5 to 
35 |a,m. As before, the volume fraction of the two constituents, bulk 
silica and nanoporous silica, are kept constant at 80% and 20%, 
respectively, for all cases considered here. 

Figure 3a and c shows the stress-strain plot comparing a two- 
hierarchy bone-like structure to a three-hierarchy self-similar assem- 
bly. We make two important observations. First, as the hierarchy 
level is increased we find a decrease in the fracture strength. 
Second, we observe an increase in defect-tolerance. Figure 3b and 
d shows the stress-strain plot comparing a two-hierarchy biocalcite- 
like structure to a three-hierarchy dissimilar structure. In this case, 
the transition from a two-hierarchy to three-hierarchy system shows 
only a small change in fracture strength, but as before an increase in 
the defect-tolerance. An important result from this analysis is a gen- 
eral trend, the observation of an increase in defect-tolerance size scale 
with the increase in number of hierarchy levels, even though the 
fracture strength tends to decrease. In Fig. 4 we take a closer look 
at the origin of the enhanced defect-tolerance. To achieve this we 
consider several samples of the dissimilar three-hierarchy structure 
with different crack sizes and identify their fracture strengths 
(Fig. 4a). We find that despite a 300% increase in crack size from 
6 |xm to 18 \xm there is only a 24% drop in fracture strength. 
A detailed analysis of the mechanics of the stress-strain curves is 



presented in Figs. 4b and 4c. We observe that the later part of the 
rising stress region before catastrophic fracture is associated with the 
opening of a large number of micro-cracks throughout the sample. 
Once these microcracks start moving and link up to the pre-crack, 
unstable crack-propagation sets in, leading to a peak in stress and 
defining the fracture strength. This effect can also be measured 
through the total new surface area created during the diffuse micro- 
cracking regime while the main crack does not move (Fig. 4c-d). 

To observe whether this pattern of increase in defect-tolerance 
continues with the number of structural hierarchies, we create a 
model with four levels of hierarchy, with results shown in Fig. 5. 
Here, the second structural level is biocalcite-like, while the third 
and fourth levels are bone-like designs. The overall volume fraction 
of the bulk-silica constituent is maintained at 80% as in all previous 
models. The overall sample size is 108 |im by 140 |im in the X-Y 
plane, and 100 nm in the Z-direction. Figure 5b shows the stress- 
strain curves for the four-hierarchy structure with various crack sizes 
from ~6 \xm to =64 |xm. Notably, almost no change in fracture 
strength is observed over this rather large change in the crack size. 
This finding confirms that the defect-tolerance has improved quite 
substantially in the four-hierarchy system over that of the two- 
hierarchy and three-hierarchy structures. 

The effect of different hierarchical structures on the fracture prop- 
erties can be captured effectively through measuring changes in the 
material crack-propagation resistance curve, the so-called R-curve. 
The R-curve measures the resistance to fracture G(Aa) for a given 
amount of crack advance length, Aa (details on the calculation of the 
fracture resistance using the J-integral are given in the Methods 
section). Figure 5c depicts R-curves for bulk silica and two, three 
and four hierarchy level structures. We observe a distinct effect of 
the addition of hierarchies on the R-curve behavior of the material. 
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Figure 4 | Detailed analysis of the source of defect tolerance in three-hierarchy structures, a. Stress-strain behavior and fracture strengths for several 
samples of the dissimilar three-hierarchy structure with different crack sizes. We find that with a 300% increase in crack size from 6 \\ra to 18 |im there is 
only a 24% drop in fracture strength, b and c. The latter part of the rising stress regime before fracture occurs involves the opening of many micro-cracks 
throughout the sample (crack size 1 1 nm). These micro-cracks are shown in red color in panel b, with the numbers indicating applied strain values. 
Once these microcracks start moving and link up to the pre-crack, there is unstable crack-propagation leading to a drop in stress and thus, the fracture 
strength. Panel d shows that this effect can also be measured through the total new surface area created during the diffuse micro-cracking regime, during 
which time the main crack remains stationary. Data shown in panels b-c for the case of a pre-crack length of 11 \\ra. 
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where both the absolute values of G versus crack advance as well as 
the slope of the R-curve, 8G/dAa, increase with the number of 
hierarchies. The slope of the R-curve has a close relation to the 
concept of defect-tolerance. As seen in Fig. 6a, in a material with a 
rising R-curve, unstable crack advance and thus catastrophic failure 
occurs when the following conditions are satisfied: 

Gapp>G/c (initiation of stable cracking), and (la) 

Gann > G and . ^^'^ > — (initiation of unstable cracking); (lb) 
" oAa oAa 

where Gapp is the applied energy release rate, G = G(Aa) is the 
material fracture toughness from the R-curve, Aa is the advance in 
the crack length and Gjc is defined as the fracture initiation tough- 
ness, where the R-curve intersects Aa = 0. This, along with the load- 
to-energy-release relation for an edge crack given by Gapp —a^na/E 
(where a and E are the applied stress and elastic modulus, respect- 
ively), implies that the load at which a certain crack size a will prop- 
agate unstably resulting in catastrophic failure can be directly 
calculated by marking off the crack size a on the negative X-axis of 
an R-curve, and constructing the tangent to the R-curve passing 
through this point, as shown in Fig. 6a. The slope of this curve is 
proportional to the load or fracture stress at which this crack propa- 
gates in an unstable fashion. (And similarly, the slope of the straight 
line that goes through Gjc is proportional to the load or fracture stress 
at which this crack starts to propagate in a stable fashion.) 
Importantly, this also implies that for a rising R-curve, the higher 
the rate of rise of G with crack advance Aa, the lesser the sensitivity of 
fracture stress with respect to crack size, and thus the higher the 
defect-tolerance length scale. Defect-tolerance length scales are thus 
closely linked to not only the absolute values of fracture-crack 



initiation values Gj/'' but also to the rising part of the R-curve, and 
directly visible through the slope of the R-curve. The key finding 
established in this paper is that these important features can be 
directly controlled by the incorporation of hierarchical structures, 
without the need to introduce new materials or any additional 
mechanisms such as interfacial processes. 

Indeed, this theoretical concept is confirmed through our own 
computational results shown in Figs. 6b and 6c. Here, the fracture 
stress is measured for two-, three- , and four-hierarchy structures for 
different crack sizes. The loss of strength as a percentage of the 
strength of the no-cracked samples has been plotted in Fig. 6b for 
all the levels of hierarchy. We observe that the sensitivity of fracture 
strength to crack size drops significantly with increasing hierarchy 
levels. Assuming that a 10% loss in fracture strength defines the 
defect-tolerance length scale. Fig. 6c confirms that this length scale 
increases non-linearly with the number of hierarchies and shows an 
strong exponential rise of the defect tolerance with increasing levels 
of hierarchies. We emphasize that this result is generic and inde- 
pendent of the specific definition of the defect-tolerance length scale. 

Discussion 

The most important result of our study is the demonstration of a 
design paradigm that explains how materials with superior mech- 
anical properties can be obtained despite the use of a single and 
structurally "poor" {i.e. brittle) base material, realized solely through 
the use of hierarchical structures without any additional mechanisms 
such as interfacial effects. Our findings may explain why biological 
materials universally show hierarchical structures, establishing that 
they provide a path to turn weakness into strength, and thereby a 
strategy for an organism's survival despite the availability of poor 
construction materials. We achieved this by examining the fracture 
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Figure 5 | Geometry of the four-hierarchy level material, stress-strain plots, and R-curve behavior, a. Geometry of two-hierarchy, three-hierarchy and 
four-hierarchy structures for comparison, with the four-hierarchy structure having a second hierarchical level that is biocalcite-like, while the third 
and fourth levels are bone-like. In the four-hierarchy structure, the color scheme is: bulk silica-red, nanoporous silica-green and blue, to show the four 
levels more clearly. The overall volume fraction of the bulk-silica constituent is kept constant at 80% in all cases, b. Stress-strain curves for the four- 
hierarchy structure with various crack sizes from =6 ^m to =64 \im. Almost no change in fracture strength is seen over this very large change in crack size 
for the four-level hierarchy material. This directly shows that the defect-tolerance has increased substantially over two-hierarchy and three-hierarchy 
structures, c, R-curve behavior of bulk silica and for two, three, and four levels of hierarchy structures (for initial crack sizes of 6.96 |im, 6.96 |im, 16.7 |.im, 
and 63.8 |^m, respectively). The R-curve measures changes in fracture toughness as a crack propagates through the change in energy released per unit 
length of stable crack advance. 
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mechanics of hierarchical structures with up to four levels of hier- 
archy, using an atomisticaUy- informed in silico mesoscale model that 
provided direct access to the fracture mechanics of hierarchical mate- 
rials. Through a series of simple computational experiments we 
demonstrated that the defect-tolerance length scale directly increases 
with number of hierarchy length scales (Figs. 3-5), quickly reaching 
—60 |xm for four levels of hierarchy (Fig. 6); presenting a significant 
departure by several orders of magnitude from the defect- tolerance 
length scale of the base material silica alone, which is on the order of a 
few nanometers. The paradigm illustrated here facilitates the construc- 
tion of materials that feature stable crack advance for a very large 
range of crack sizes, and without catastrophic material breakdown. 

Notably, stable crack advance is not catastrophic or disabling for 
biological materials, since materials such as bone have the ability to 
self-heal over time from fracture^' or from fatigue-induced 
damage*"'". This agrees directly with our finding that hierarchical 
materials feature a rising R-curve behavior, across several micro- 
meters and larger length scales, which promotes stable crack advance 
in the material even if it is subjected to large loads and allows it to be 
useful beyond its fracture initiation load point. Large loads are dis- 
sipated in multiple-hierarchy materials through the initiation and 
arrest of cracks that occur at different length scales (see Figs. 2d, 
4b and 4c). This is in contrast to single hierarchy materials that would 
shatter with a single crack propagating through the material. If the 
load does not increase to levels that would lead to unstable crack 
propagation, this mechanism provides a tissue or an organism the 
necessary time to repair the material in order to restore its full load 



carrying capacity. Our data confirms that a higher number of hier- 
archies provide an improved rising R-curve behavior despite no 
additional material used or new mechanisms introduced. 

Future work could be aimed at optimizing the R-curve behavior 
over several micrometers or higher length-scale of crack advance. 
Both the number of hierarchies, and the design of the geometry of 
individual hierarchy levels, would have different effects on the entire 
R-curve shape. Our mesoscale model presents an effective tool that 
can be utilized for systematic design optimization to maximize 
R-curve magnitudes and R-curve slope improvements. The model 
could also be extended to larger length scales, potentially sub-mm, 
such that R-curve toughness improvements can be seen over length 
scales that approach macroscopic sample sizes. Further work could 
also be focused on an improved modeling of the interfaces between 
soft/hard phases at each hierarchy level, for example to include a 
specific constitutive behavior. This would enable one to capture 
the contributions of the interfaces to the toughness at various levels 
and open additional avenues to improve the material performance. 
However, it is emphasized here that the model used in this article was 
chosen deliberately to preclude any additional mechanisms at inter- 
faces to provide a clean model that examines effects solely due to 
hierarchical structures. Despite the mechanistic insight developed 
here, the model used in this study is qualitative in nature and the 
toughness values predicted have to be carefully validated against 
experimental data. Nevertheless, our results show a manifold 
increase of key fracture parameters which should hold regardless 
of the details of the model used. 
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Figure 6 | Relation between the R-curve, rate of fracture strength change with crack size, and the defect-tolerance length scale for varied number of 
hierarchies, a, Definition of variables a (initial crack length) and Aa (crack advance length), b. Schematic of R-curve G(Afl) in a material with a rising 
R-curve resistance and Unk to unstable crack propagation. The rising R-curve shown here resembles those seen in hierarchical structures (e.g.: Fig. 5a). 
The load at which a crack with a certain size will propagate unstably and thus cause catastrophic failure can be calculated by marking off the crack size on 
the negative-X axis and constructing the tangent to the R-curve passing through this point. The slope of this curve is proportional to the load or fracture 
stress at which this crack propagates unstably. For comparison the plot also shows a second R-curve (G'(Afl)) that does not rise as rapidly, resembling 
those R-curves seen in Fig. 5a for fewer levels of hierarchies. The fracture stress at which the crack propagates unstably is much smaller, indicating a smaller 
defect tolerance, b. Fracture stress as a percentage loss from the strength of structures with no cracks, measured for two-, three- , and four-hierarchy 
structures and for varied crack sizes. The shaded region shows the crack sizes with less than 10% loss in strength. The data shows that the sensitivity of the 
fracture strength to crack size decreases substantially with increasing hierarchy level, c. Plot of the defect-tolerant length scale over the number of 
hierarchies. The defect-tolerant length-scale reaches —60 [im with four levels of hierarchy. The red line represents an exponential fit of I = 1.17 
exp(1.255N), where L is the defect-tolerant length scale in micrometers, and Nis the number of hierarchy levels. The limiting value of one level of 
hierarchy, i.e., bulk silica, features a vanishing defect-tolerant response on the order of nanometers. 
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Methods 

We use a multiscale bottom-up computational methodology to study the effect of 
hierarchical design on material properties. At the micron length scale, we use a simple 
mesoscopic model that describes the material as a collection of particles connected by 
complex nonlinear springs^^, with all parameters informed based on full- atomistic 
simulations based on reactive force field modeling. The model is designed by our 
desire to develop a simple model to derive generic insight into the mechanical 
properties and deformation mechanisms valid for a broader class of hierarchical 
biological materials. We note that even though such a simple model formulation does 
not allow us to derive quantitative conclusions for phenomena pertaining to specific 
types of biological materials, it enables us to understand universal, generic relation- 
ships between underlying molecular mechanisms, resulting nonlinear properties of 
the material, and the failure behavior of hierarchical materials. The level of abstrac- 
tion used in the present study is critical for our goal to establish fundamental insight 
into a broader class of materials. 

At the nanoscale we use full- atomistic molecular dynamics simulations to study the 
mechanics of the nanoporous silica structures. We use the first principles based 
reactive ReaxFF atomistic force field^^'*^, which effectively captures the energy land- 
scape associated with chemical bond changes, including bond breaking. The ReaxFF 
description is based on a bond-length bond-order description and fitted to density- 
functional calculations of energy landscapes of bond-distortion, breaking and 
forming events. A variety of Si-O clusters were used for fitting parameters, as also 
energetics of bulk crystalline phases of silicon and silica under tension and com- 
pression; all obtained from quantum mechanics calculations (for details see*^'*^). The 
ReaxFF potential has been used successfully in predicting fracture phenomena in 
silicon and silica, and has been validated against experimental data^^'^-*^. At the 
micron length scale, with the mesoscale spring-lattice network model, structural 
loading is carried out by stepped edge displacement boundary conditions and 
relaxing the global positions of all material particles using a conjugate gradient energy 
minimization scheme*^. 

The fracture toughness in the mesoscale model is calculated by invoking the 
J-integral^^. R-curves are plotted by calculating structure fracture toughness for dif- 
ferent values of stable crack advance^^. Stable crack advance for every load config- 
uration is noted by fmding the crack tip location. A particular bond is regarded as 
broken when its deformation exceeds the cutoff for the interaction. Crack surfaces are 
visualized by finding all spring bonds which have snapped for a given load. The 
J-integral is used to fmd the energy release rate for a given amount of stable crack 
advance. The plot of the J-integral from the start of crack initiation throughout crack 
propagation provides the R-curve*^ for the material, i.e. how its fracture toughness 
changes as a function of stable crack advance. J-integral calculations for different 
structures are done with large initial crack sizes. In all cases the initial cracks are 
chosen large enough to prevent a diffuse micro -cracking response in the entire 
material, as this would prevent us from carrying out the calculation of the J-integral 
through the method demonstrated above, due to unavailability of an integration 
region free of cracks. 

A more detailed description of the computational methods is provided in the 
Supplementary Information section. 
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